Discontinuous Transition of a Multistage Independent Cascade Model on Networks 
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We study a multistage independent cascade (MIC) model in complex networks. This model is 
parameterized by two probabilities: T\ is the probability that a node adopting a fad increases the 
awareness of a neighboring susceptible node until it abandons the fad, and Tb is the probability 
that an adopter directly causes a susceptible node to adopt the fad. We formulate a framework of 
tree approximation for the MIC model on an uncorrelated network with an arbitrary given degree 
distribution {pk}- As an application, we study this model on a random regular network with 
degree k — 6 to show that it has a rich phase diagram including continuous and discontinuous 
transition lines for the percolation of fads as well as a continuous transition line for the percolation 
of susceptible nodes. In particular, the percolation transition of fads is discontinuous (continuous) 
when Ti is larger (smaller) than a certain value. Furthermore, the phase boundaries drastically 
change by assigning a finite fraction of initial adopters. We discuss this sensitivity by both tree 
approximation and numerical simulation. 

PACS numbers: 89.75.Hc,87.23.Ge,05.70.Fh,64.60.aq 



INTRODUCTION 

Word-of-mouth phenomena often catapult books, 
movies, and music out of obscurity and into popular- 
ity. Similarly, we witness similar propagative effects in 
voting behavior, the spread of rumors, and diffusion of 
innovations. To explain the mechanism causing such in- 
formation propagation, Watts proposed a network model 
(the Watts model) in which the decision made by each 
node (person) to adopt an innovation is determined by 
the relative number of its adjacent nodes (friends) who 
have already adopted it In this model, the system 
transits from the global cascade phase, where some small 
initialized shocks are propagated though a whole system, 
to the local cascade phase, when the threshold of their 
decision making increases. This transition is discontin- 
uous: the final fraction of adopters rises discontinuously 
and its distribution over many trials is bimodal even at 
the transition threshold. Following the seminal paper by 
Watts, several authors have investigated the effects of 
network topology on the transition threshold 0-0] • 

The Watts model assumes that a node is permanently 
active in the sense that the node will continue to adopt 
an innovation once it has done so. On the other hand, 
the susceptible-infected-removed (SIR) model @ and the 
independent cascade (IC) model [§] (including the SIR 
model with transmissibility fiol ]) have been investigated 
as models of transient fads. In the IC model, an adopter 
of a fad makes only one attempt, with a certain probabil- 
ity, to persuade each susceptible neighbor to adopt it, and 
after that attempt, the adopter abandons the fad. The 
SIR model and the IC model have been extensively stud- 
ied both analytically and numerically on complex net- 
works [Icil - fl2| . On such networks, these models do not 
show a discontinuous transition, but a continuous one 



In recent years, several multistage contagion models 
have been reported Even when a fad is transient, 

it may show a discontinuous transition when the number 
of possible states for each node is increased. Krapivsky 
et al. proposed a model of transient fads (the fad model) 
by incorporating social reinforcement into the SIR model 
jig . Social reinforcement means that each node requires 
multiple prompts from adjacent adopters before adopt- 
ing a fad. In this model, a susceptible node takes one of 
M awareness levels 0, 1, • • • , and M — 1, and it becomes 
an adopter when the level reaches M. Each adopter in- 
creases the awareness level of the adjacent neighbors by 
one at unit rate until the adopters begin to abandon 
the fad at a certain rate. Note that the case of no re- 
inforcement (M = 1) coincides with the standard SIR 
model. Krapivsky et al. studied this model in a mean- 
field framework to show that the transition is continuous 
for M = 1 and discontinuous for M > 2. However, noth- 
ing is known about how this model behaves when placed 
on a network. 

In this paper, we propose a simple multistage indepen- 
dent cascade (MIC) model, inspired by the fad model. 
This model is parameterized by two probabilities: T± is 
the probability that a node adopting a fad increases the 
awareness level of a neighboring susceptible node until 
the node abandons the fad, and T2 is the probability in 
which an adopter directly persuades a susceptible node 
to adopt the fad. We formulate a framework of tree ap- 
proximation for the MIC model on an uncorrelated net- 
work with an arbitrary given degree distribution {pk}- 
As an application, we study this model on a random reg- 
ular network with degree k = 6 to show that it has a 
rich phase diagram. The percolation transition of fads 
is discontinuous (continuous) when T± is larger (smaller) 
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FIG. 1: Transition rules of the MIC model. The values on 
the lines indicate the transmission probability. The events 
indicated by the dashed lines (0 — ¥ 1, — ¥ 2 and 1 — > 2) occur 
when the node is adjacent to an adopter, whereas the event 
indicated by the solid line (2 — ¥ 3) occurs irrespective of the 
neighbor's state. 

than a certain value. Furthermore, the phase boundaries 
drastically change by setting a finite fraction of initial 
adopters. We discuss this sensitivity by both tree ap- 
proximation and numerical simulation. 



PHASE DIAGRAM OF THE MIC MODEL IN 
THE p -> LIMIT 

In this section, we consider the dynamics starting from 
an infinitesimal fraction of seeds (p = 0+) on uncor- 
rected networks with a given degree distribution {pk}- 
Newman analyzed the SIR model with a transmissibility 
T (the MIC model with T t = and T 2 = T) by using 
generating functions [HI EH- According to his result, the 
system experiences two phase transitions in increasing T: 
the giant component of abandoners appears at the per- 
colation threshold of the S3-network T c = (k) / (k 2 — k) 
((•) means the average of a quantity weighted by pk), and 
the giant component of susceptible nodes disappears at 
the percolation threshold of the So-network T c i, Corre- 
spondingly, there exist three distinct phases: S™ a * > 
and S^ ax = below T c , S^ x > and S;f ax > for 
T c <T <T C >, and S nax = and S 3 nax > above T c , . 
In the following, we generalize his method to the MIC 
model. 



MODEL 

Let us consider a network with N nodes. Each node 
takes one of the four states: (0) susceptible (without 
awareness), (1) informed (susceptible with awareness), 
(2) adopter of the fad, and (3) abandoner of the fad. 
The dynamics of this model is given as follows (Fig|T]): 

(i) Select randomly the fraction p of nodes as the initial 
adopters (seeds). Set the state of other nodes as 
being susceptible. 

(ii) For each adopter, change each of the neighboring 
susceptible nodes to be cither informed (0 — >• 1) 
with probability T\ or an adopter (0 — > 2) with 
probability T 2 and change each of the neighboring 
informed nodes to an adopter (1 — > 2) with proba- 
bility T = Ti + T 2 (< 1). 

(hi) Change all the old adopters to abandoners (2 — > 3). 

(iv) Repeat (ii) and (hi) until no adopter exists in the 
network. 

Note that the case with T\ = coincides with the stan- 
dard IC model (the SIR model with transmissibility T2), 
and the case with T 2 — is regarded as an IC version of 
the fad model (27j|. 

In the following sections, we denote the fraction of 
nodes in states 0, 1, 2, and 3 in the final state by So, 
Si, S2, and 53, respectively. Also, we designate the net- 
works composed of susceptible nodes and abandoners the 
So-network and the ^-network, respectively. The frac- 
tions of the largest components of the So-network, and 
the 53-network are denoted by S™ 11 * and S™ a *, respec- 
tively. 



Percolation of the ^-network 

Let u be the probability of a susceptible node not being 
informed by the prompting of a randomly chosen neigh- 
boring (RCN) node, and v be the probability that a RCN 
node does not adopt the fad in the final state of dynam- 
ics. Then u and v are related as 

u = v + (l-T)(l-v) = l-T + Tv. (1) 

To evaluate these values, it is convenient to introduce the 
generating function of the degree distribution {pk} as 

F (x) = 5> fc x fc , (2) 

k 

Fi( X ) = ||f| = ^E(^lW^- (3) 

By using these functions we can express the probabilities 
of a node being susceptible and informed, which are equal 
to So and Si respectively, as 

So=F (u), (4) 
Si=Ti(l-v)F^u), (5) 
S 3 = 1 - So - Si- (6) 

The last equation follows from the fact that S 2 = at the 
final state. Now let us divide v into two parts, v = vo+vi, 
where vo and vi are the probabilities of a RCN node 
being susceptible and informed, respectively. The self- 
consistent equations for vo and vi are then given as 

v Q = Fi(u), (7) 
vi=Ti{l-v)F[{u), (8) 
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or by adding them we have 

v = F 1 (u)+T 1 (l-v)Fl(u), (9) 

which can be expressed in terms of u as 

u = l-T + TF 1 {u) + T 1 {l-u)F[{u). (10) 

Let us analyze in some detail the solution of this equation 
in the following. 

For any values of T\ and T 2 , there exists a trivial so- 
lution u = 1, meaning So = 1, S\ = 0, and thus S3 = 0. 
To find nontrivial solutions we solve ([TO]) with respect to 
T 2 as 



x , A + TiFUu) 

r 2 = yo (»;ri) = (i-«) -T t . 



(ii) 



This gives T 2 such that a given value of u < 1 is a solution 
of (fini) with a fixed T x . For small u = 1 - u < 1, (HH) is 
expanded to give 

1 caO-ciTi) 

T 2 = 1 r-a U 

ci 2cf 



(3c 2 - 2cic 3 ) + ci(4cic 3 - 3cg)gt _ 2 
12c? U 



where c„ = i<i (1). Taking the limit u — > 1— (u —> 0+), 
we have 



(12) 



T 2 =T C = cA = 



1 



(k) 



F{(1) (fc 2 -fc)' 



(13) 



above which the trivial solution u = 1 becomes unstable 
and S3 should be nonzero. 

For Ti < c^ 1 = T c fixed, a nontrivial solution exists 
only for T 2 > T c and S 3 -> when T 2 -> T c +, so that 
T 2 = T c is the continuous transition point of S3. For 
T\ > T c , on the other hand, u < 1 is possible even for 
T 2 < T c . To see this let us rewrite (|T2J| for small A x = 
T\ — T c and A 2 = T c — T 2 , which has the structure as 



A 2 = 2aAiw - bu 2 



(14) 



with some positive constants a and b. For Ai > 0, A 2 
can be positive for small u > 0, and thus there are two 
branches of the nontrivial solution for Tj 1 < T 2 < T c with 



J- ,1 



T - —A 2 



(15) 



and the lower branch vanishes at T 2 = T c while the upper 
branch survives above T c giving a positive S3. For T 2 < 
Tj" , on the other hand, a nontrivial solution does not 
exist and the only possibility is S3 = 0. This means 
that S3 should jump dis continuously from zero to the 
upper branch at somewhere < T 2 < = T c . In 
general, depends on T\ and is given by T^(Ti) = 
min <£ (u;Ti). 

0<ti<l 



Percolation of S'o-network 

When T 2 increases with T± fixed, the So-network grad- 
ually shrinks and finally breaks into a lot of finitely con- 
nected components at T 2 = T 2c ' ■ Here, let us examine 
this percolation problem. 

The probabilities that a RCN node of a susceptible 
node in the So-network are susceptible, informed, and an 
adopter are Vq, V\, and (1— T)(l— v), respectively. There- 
fore, the generating function for the degree distribution 
of the So-network is given as 



G (x) 



F (v x + Vl + {I - T)(l - v)) 
Fo(u) 

F (v (x - 1) + u) 
F (u) 



(16) 



Similarly, the generating function for the excess degree 
distribution of the So-network is derived as 

r M _ G' (x) _ F^vojx-^+u) 



G' (l) 



Fi(u) 



Because the mean excess degree of the So-network is 
given as G^l), the value of u at the percolation threshold 
(Ti c /,T 2c /) should satisfy 



l = G' 1 (l) = v 



F{(u) 
Fi(uY 



(18) 



In the percolating phase the largest component S™ 8 * is 
given by 



S max = So (1 - G (vs)) . 



(19) 



where vs is the probability that an So-cluster connected 
via a randomly chosen link to a susceptible node is finite, 

vs = G 1 (v s ). (20) 



Phase diagram for random regular network 

As an application, we consider the phase diagram of 
the MIC model on a random regular network with degree 
k = 6. By applying the present tree approximation, we 
have the phase diagram in T\-T 2 space as shown in FigJ5] 
The transition line of the S3-pcrcolation is independent of 
Ti and located at T 2 = T c = l/(fc-l) = 0.2 (the vertical 
line). However, the order of transition is continuous if 
T\ < T c , and otherwise discontinuous. In Fig|3l we plot 
S3 with several values of T\ as a function of T 2 . For large 
Ti, S 3 is multivalued in T d L < T 2 < T% = T 2c . Since 
dynamics starts from an infinitesimal fraction of seeds, S3 
remains as the lower branch below Tj 7 and jumps to the 
upper branch at T 2 = , implying that a discontinuous 
occurrence of the giant component of the S3-network at 
T 2 = Tj 7 . Such a jump disappears when T\ < T± c . 
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FIG. 2: Phase diagram of the MIC model with p = 04- . The 
thick-solid and thick-dashed lines denote the continuous tran- 
sition line for percolation of the S3-network T^d and of the 
So-network Tbc, respectively. The thin-solid and dashed lines 
denote the upper (Tj 7 ) and lower (T<f) bounds, respectively, 
between which (|10|) has two nontrivial solutions. 



FIG. 4: Phase diagram of the MIC model with p = 0.01. 
The thin-solid and dashed lines denote the upper (Tj 7 ) and 
lower (Ta) bounds, respectively, between which (|10p has two 
nontrivial solutions. The thick-dashed line denotes the tran- 
sition line of percolation of the So-network T 2c < . The marks 
denote the transition point of percolation of the S^-network 
obtained by Monte-Carlo simulations. 




FIG. 3: Analytical result for S 3 with Ti = 0, 0.2, and 0.4 
(from right to left) as a function of Tb. T% and indicate 
the lower and upper bounds, respectively, between which (|10p 
has two nontrivial solutions. 



The region of coexistence of the giant components of 
the Sq~ and ^-networks exists in T 2c < T% < Tic'iTi) 
for small T\ . This region shrinks with increasing T\ , and 
T2c'(Ti) reaches T 2c at a certain value of T\ above which 
the percolation of the S^-network and the disappearance 
of the giant component of the So-network occur simulta- 
neously at T 2 = T 2c . 

The tricritical point (Xi c ,T2 C ) of the random regular 
network with degree k is given by (l/(fc — 1), l/(fc — 1)). 
We note that {Ti c ,T 2c ) = (0.5,0.5) when k = 3, so there 
is no room for T\ > T& at T 2 = T 2c since T = T\ + 
T2 < 1. This shows that the MIC model on the random 
regular network with degree k = 3 does not show the 
discontinuous transition of the 53-network, but only the 
continuous one. 



SENSITIVITY TO THE FRACTION OF INITIAL 
ADOPTERS 

Tree approximation 

On the basis of the formalism in the previous section, 
it is easy to find the self-consistent equations for the case 
where the fraction of initial adopters is p > 0. Let p = 
1 — p be the initial fraction of susceptible nodes. Then 
for the present model v = vq + V\ where 

v = pF^u), vi = pT x (l - v)F{(u), (21) 

the relation u = 1 — T + Tv is unchanged. In terms of u, 
we thus obtain 

u = 1 - T(l - pF^u)) + pT^l - u)F[{u). (22) 



By using the solution, we obtain 



So^pFaiu), Si=p 



Ti(l 



T 



(23) 



and S3 — 1 — 5b — Si . 

For p > 0, percolation transition of the ^-network, 
from zero to nonzero S l 3 Tlax , can take place as long as p 
is smaller than the percolation density (see the next sub- 
section). In contrast, there is no trivial solution of (|22|) 
for all regions of parameters and no continuous transi- 
tion from zero to nonzero S3 because all seeds will be- 
come abandoners, so S3 > p. Nevertheless, discontinuous 
transition of S3 is still possible. Equation (ITU is modified 
as 



T 2 = ip p {u;Ti) = (1 - u) 



l + p^Fjiu) 
l-pF 1 {u) 



(24) 
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FIG. 5: Comparison of S3 with the numerical result for 
T\ — 0.4 and p = 0.01. The analytic result (solid curve) 
and the result obtained by simulation with TV = 1000(x), 
2000(+), 4000(D), 8000(H), and 16000(0) are shown. Here 
£3/(1 — S3) is plotted instead of S3 itself to see the behaviors 
near S3 = and 1. The analytic evaluation gives T% ~ 0.1306 
and T% ~ 0.1587. 



FIG. 7: Monte Carlo results for S^ with Ti = 0.25 as a 
function of T 2 . Data for N = 1000(x), 2000(+), 4000(D), 
8000(B), and 16000(0) are shown. The solid line is drawn 
by CES), and gives T 2c , ~ 0.22576. 



Results for random regular network 
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FIG. 6: 

p = 0.01 as a function of T 2 . Data for TV = 1000(x), 2000(+), 
4000(D), 8000(B), and 16000(0) are shown. The discontin- 
uous transition point of S3 and the transition point of S3- 
percolation are estimated from these data as T2 ~ 0.158 and 
T2 ~ 0.106, respectively. 



and the lower bound and the upper bound T% of the 
discontinuous transition point can be given by the local 
minimum and local maximum of (p p (u;Ti) in < u < 1, 
respectively. We note that there is a value of p, depending 
on Ti, above which <p p (u; Ti) is monotonically decreasing 
and these bounds do not exist (see Fig[5]). 

For the So-network, the formulation for deriving the 
percolation threshold T c > and the ratio of the largest com- 
ponent S™ 3 * is unchanged, so that (JTSJ) and (JTHJl still hold 
with (dU) and (|2"5]), 



Now let us discuss the sensitivity of the fraction of 
the initial adopters by considering the random regular 
network with degree 6. Here, we performed Monte Carlo 
simulations of the MIC model to verify our analytical 
results. The number of nodes used is N = 1000, 2000, 
4000, 8000, and 16000. The number of graph realizations 
is 100, and the number of trials on each graph is 100. 

In FigpH we show the phase diagram for the case with 
p = 0.01. The upper and lower bounds of the discon- 
tinuous transition of the S^-network and the percolation 
threshold of the So-network are drawn by the tree ap- 
proximation. We find that the upper bound shifts to the 
left from the vertical line (the discontinuous transition 
line of the case with p = 0+). Wc also plot numer- 
ical results for the percolation transition points of the 
S^-network (our tree approximation does not give S'™ ax 
when p > 0). Each transition point is estimated from 
the cross point of the fractal exponent ljl, 2(J , which is 
defined as 



Urn lo gAr ST x (iV) 



A' 

log^ 



(2N) -logS^ ax {N/2) 
log 2N - log N/2 



max 
3 



(25) 



Note that the crossing point of the fractal exponent indi- 
cates that the transition is continuous. We find that the 
continuous transition points of the S3 network are outside 
the interval between the lower and upper bounds. This 
means that the percolation of the S^-network and the dis- 
continuous jump of S3 occur at different points: when T2 
increases with T\ fixed, the continuous transition of the 
^-network occurs first, and then a discontinuous change 
in S3 occurs at a larger value of Ti between the lower 
and upper bounds. 



6 




FIG. 8: T 2 dependence of S 3 with Ti 
and 0.1 (from left to right). 



0.4 and p = 0, 0.01, 



the phase diagram of the MIC model on a random reg- 
ular graph and confirmed that the results support well 
our analytical predictions. 

Our tree approximation for the present MIC model is 
valid only for uncorrclated networks, and it is beyond our 
scope to investigate how the cooperative behavior of this 
model depends on the properties of the network on which 
the dynamics is performed (e.g., clustering, degree cor- 
relation, and community structure). For example, clus- 



tering decreases the threshold of the SIR model |21| , but 



an innovation in the Watts model is easier to propagate 
in rather clustered networks (as compared to random 
networks) @, |(|. How does the MIC model behave in 
clustered networks? In order to answer such questions, 
further research is needed. 



In FigjSl we compare the analytical result and nu- 
merical evaluations for S3 with T\ = 0.4. We find 
that the analytical result coincides well with the numer- 
ical one except near the discontinuous transition point 
T 2 = Tj^Ti) ~ 0.1587. Near the transition point, nu- 
merically evaluated S3 approaches the analytical curve 
when N increases. Although it is not analytically con- 
firmed that the upper bound Tj 7 coincides with the ac- 
tual discontinuous transition point when the fraction of 
the initial adopters is finite, the jump of S3 from the lower 
branch to the upper branch occurs very near or just at 

1 d ■ 

We plot the Monte Carlo result for S^ with Ti = 0.4 
and S£ lax with Ti = 0.25 for p = 0.01 in Figs© and [0 
respectively. From Fig|S]we find that the S^-network first 
percolates but grows discontinuously at a larger value of 
T-2- In Fig 13 we have a good correspondence between the 
numerical results for S™ ax and Eq.(fl"9|). 

Finally, we show analytical curves of S3 with several 
values of p and a fixed T\ = 0.4 in Figj8] When p is small, 
we observe the precursor of the discontinuous transition 
of S3 in the same way as in Fig|3] However, such a be- 
havior disappears with increasing p. The system does 
not show the discontinuous transition if the fraction of 
the initial adopters p is sufficiently large. 



SUMMARY 

We have investigated a simple multistage independent 
cascade (MIC) model on networks. Under the local tree 
approximation, we have derived the formulation for eval- 
uating Sq, Si, and S3 and thus the phase diagram. In 
particular, percolation of the adopters occurs discontin- 
uously for large values of T\. Moreover, when dynamics 
starts from a finite fraction of initial adopters, the discon- 
tinuous jump in the number of adopters occurs after the 
adopters percolate. As an application, we have studied 
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